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I. INTRODUCTION 



In heavy ion collisions at very high energies, the matter created differs qualitatively from what is traditionally studied 
in nuclear and elementary particle physics. In the initial stages of the collision, copious production of gluons and 
quarks in a large volume leads to a rapid increase in the entropy, and the distinct possibility of a new phase of 
matter characterized by deconfined degrees of freedom. One therefore hopes that relativistic heavy ion experiments 
can provide insight into the structure of the QCD vacuum, deconfinement, and chiral symmetry restoration. 

The hot transient system eventually evolves into a gas of hadrons at high energy densities, whose properties may 
be studied theoretically using, for example, hadronic cascades [|j]-f3j. In principle, these models provide information 
on the early, dense phase by tracing the evolution of the system from hadronization to freeze-out. Of course, in 
ultrarelativistic heavy ion collisions, most the of the produced secondaries are pions. For example, in central Au+Au 
collisions at center-of-mass energies of 200 A GeV estimates from the FRITIOF event generator suggest that ~ 4000 
pions per isospin state might be produced. Further, recent measurements [Q at lower energies and comparison to 
simulations || show that freeze-out source sizes probably deviate quite drastically from a simple multiplicity scaling 
law: present calculations indicate 10-20 fm Au+Au source radii at y/s = 200 A ■ GeV. In any event, these high 
energy collisions might well create highly degenerate bose systems, and even possibly Bose-Einstein condensates 
(BEC) . Since practical conclusions from dynamical simulations |(| depend qualitatively on the effect of the medium 
on particle interactions |?],|| , one needs to better understand the properties of such degenerate systems of pions within 
the environment of a relativistic heavy ion collision. 

Non-relativistically, the problem of interacting, degenerate bose systems has been discussed extensively by several 
authors. Evans and Imry |J established the pairing theory of a bose superfluid in analogy to the BCS theory of 
superconductivity. For an attractive interaction, the resulting gap equation may have a non-trivial solution. Further, 
though, there appears the possibility of having a macroscopic occupation of the k — particle state when the 
corresponding BCS quasiparticle energy vanishes. In turn, this leads to a spectrum which is linear and gapless in 
the long wavelength limit |J. In a second paper, Evans and Rashid rederived the equations of Ref. |i| using 
the Hartree-Fock-Gorkov decoupling method, and solved them for the case of superfluid helium. This boson pairing 



theory has been generalized by Dorre et al. 11 1, who carried out a thermodynamic variational calculation with a 
trial Hamiltonian containing a c-number part. An extensive discussion on the boson pairing problem is also given by 
Nozieres and Saint James [[l2| . 

It has further been shown by Stoof jl3| and, independently, Chanfray et al. |l4| that the critical temperature T c for 
the transition from the normal phase to the phase with a non-vanishing gap (the Evans-Rashid transition) is given by 
a "Thouless criterion" |l5|,|l6| for the bosonic T-matrix in the quasiparticle approximation, in analogy to the fermion 
case. Moreover, it has been demonstrated that there exists a second critical temperature Tbec < T c , where the 
condition for the macroscopic occupation of the zero momentum mode of Ref. M is fulfilled ju|,[u] . The mechanism 
is the same as for Bose-Einstein condensation in the ideal bose gas jl3| . 

Here we wish to consider n-n interactions in the presence of a dense and hot pion gas along the lines of a previous 
approach |fjj|l4j]. We address the question of pion pair formation and the pion dispersion relation in a thermal 
medium, first in a qualitative way (section II), then in a more detailed numerical calculation with a realistic two pion 
interaction (section III). As we shall see in section IV, the in-medium 27r propagator exhibits a pole above a certain 
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critical temperature, signaling a possible instability with respect to pion pair formation. We conclude in section V 
with a discussion of the effect in high energy heavy ion collisions. 

The effects we present here require rather large phase space densities for the pions, but are independent of whether full 
thermal equilibration has been reached. Nonetheless, we choose to couch the discussion in terms of thermal language, 
because it is convenient, but also because the actual situation is probably not too far removed from it. Dynamical 
calculations |||| show that a high degree of thermal equilibration is quite reasonable. Chemical equilibration, on 
the other hand, may well cease at later stages of the system's evolution and lead to a condensation of pions in the 
adiabatic limit. Of course, the system actually expands rather rapidly, but nonetheless large chemical potentials 
(fi ~ 130 MeV) may be built up by freezeout (T ~ 100 MeV). One might thus expect large phase space occupation 
numbers at low momenta, which drive the pion pair formation that we discuss here. 



In order to treat the gas of interacting pions we will use the boson pairing theory of Evans et al. |5|,[l0| • In analogy to 
the fcrmion (BCS) case, they obtain a system of coupled equations for the gap energy and the density by linearizing 
the equations of motion. The usual Thouless criterion for fermions can be established analogously for the bose system, 
and yields the critical temperature below which the gap equation begins to exhibit non-trivial solutions. However, in 
contrast to the fermion case, a second, lower, critical temperature appears at which the quasiparticle energy vanishes 
at zero momentum. This temperature is associated with the Bose-Einstein condensation (BEC) of single bosons, in 
analogy to the ideal bose gas, as discussed in Ref. |l4| , and in detail for atomic cesium in Ref. |13| . An interesting 
feature of the formalism developed by Evans et al. |)) is that below the second critical temperature the dispersion 
relation for the single bosons is of the Bogoliubov form, i.e., linear, or phonon like, for small momenta. 

In this section, we illustrate these remarks concerning the Evans-Rashid transition for a pion gas in a qualitative way, 
returning to a more detailed numerical calculation in section III. While relativistic kinematics is taken into account, 
corrections from backward diagrams are ignored. We shall see in section III that such an approximation is justified for 
the physical regions in which a solution to the gap equation exists. For clarity in this preliminary discussion, we shall 
also generally neglect the fc-dependence of the self-energy, S(fc), and further condense it into the chemical potential. 

The gap equation for the pion pairs is derived in the appendix, using Gorkov decoupling: 



with e(k) = ui{k) — /z, and where uj{k) is the free pion dispersion. The coth-factor represents the medium effect for 
a thermalized pion gas at temperature T and chemical potential /i, and V(k,k',E) is the as yet unspecified bare 
two-particle interaction. The corresponding pion density is 



II. THE EVANS-RASHID TRANSITION IN A HOT PION GAS 




(1) 



where the quasiparticle energy is given by 



E(k) = vW~W 



(2) 




(3) 



3 



In spite of the formal similarities of Eq. (|l]) with the corresponding fermionic gap equation, there are important 
differences: For bosons, the A 2 is subtracted in E(k) (fermions: added), and the temperature factor is a hyperbolic 
cotangent (fermions: tanh). 

We discuss the solution to the gap equation for decreasing temperature, at a fixed value of the chemical potential a. 
The possibility of a finite chemical potential in a pion gas has been pointed out in the introduction. At very high 
temperatures, the gap equation (|l|) has only the trivial solution A = 0, and Eq. (||) is the usual quasiparticle density. 
The dispersion relation is also of the usual form 



E(k) = e(k) = uj(k) - (j, = \Jk 2 + ml - fx . (4) 

With decreasing temperature, however, a critical temperature T" is reached, at which the gap equation first 
exhibits a non-trivial solution A / 0. The value of T" may be found by linearizing the gap equation, i.e., setting 
E{k) w e(fc) in Eq. ([!]). We return to this point in section III, showing that the resulting equation for T" is identical 
to the condition for a pole in the two-pion T-matrix at the particular energy E = 2u (for total momentum K = 0). 
Thus we have a bosonic version of the well-known Thouless criterion for the onset of superfluidity in fermion systems 
with attractive interactions. 

Below the critical temperature T™ the order parameter A becomes finite, and the corresponding dispersion relation 
is now given by Eq. (||). As the temperature drops further, |A| increases to a point where the condition |A(fc = 0)| = 
\m^ — fi\ is reached. This is the maximum possible value of |A|, since otherwise imaginary quasiparticle energies 
result. It defines a second critical temperature Tj, below which the occupation no of the zero momentum state 
becomes macroscopically large because E(k) — >0 for k— >0 |J. The possibility of a macroscopic occupation of the 
k = mode below follows from the pion density Eq. (H|): for E(k = 0) = 0, the k = contribution to the density 
must be treated separately, as in the case of the ideal bose gas. A similar comment applies to Eq. (]]]) for the gap, so 
that we obtain the two inhomogeneous equations 

A(fc) = -iv(M,2A0no - \ _E V ^ k '^^r ) coih ^W L ' (5) 

fe'^O 



n + 



2E(k') 2T 2 



(6) 



In contrast to the ideal bose case, the condensation of quasiparticles happens at /i < , because of the finite value 
of the gap. Below Tjf the dispersion relation is given by 



E(k) = \/uj(k) 2 - 2w(fc)u + 2/xm^ - ml , 

« W2(m.-/i)^- + ^(-^-) 2 , (7) 

in the small k, non-relativistic approximation. Thus, instead of the usual fc 2 -behavior, the pion dispersion is linear in 
the long wavelength limit. Eq. (Q) may be rewritten in the more usual form of the well-known Bogoliubov dispersion 
relation Ca] for a weakly interacting bose gas: 



E(k) = W 1^(0,0,2^)1 n ^- + 0(fc*) . (8) 



Here, we have used — fi = — V(0, 0, 2fi) no/2, which follows from Eq. (^|) for sufficiently low temperatures. 
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III. NUMERICAL RESULTS FOR THE GAP EQUATION 



We now consider the qualitative discussion of the previous section in more detail, by numerically solving our system 
of equations for a realistic pion-pion interaction in the i = I = channel. We choose a rank-2 separable tt-tt 
interaction inspired by the linear a-model (see appendix) which possesses all the desired low energy chiral properties, 
as is explicitly discussed in Ref. [i7fl . For vanishing incoming total momentum, K = 0, it reads (see Eq. ( A.23| )) 



(k,-k\ Vr = o(£) | k',-k') 



v(k) Ml - ml 



2u)(k) 
1 

2cj(k) 



f2 

■I TV 



E 2 - mi 



4u(k)u(k') - 2ml 



E 2 



(k | M B {E) I k') 



-Ml 
1 

2w(fc') 



Ml 



v(k') 
2u{k') 



where, for later convenience, we have introduced the bare invariant matrix 

(k | M B (E) I k') = X 1 (E) Vl (k) Vl (k') + \ 2 v 2 (k)v 2 {k') 



with notation ui(fc) = v{k) 
Xi(E) = 



[1 + (k/Smn) 2 ]' 1 , v 2 (k) = (w(k)/m^)v{k), and 



Ml 



f 2 

J TV 



E 2 



2ml 



E 2 - Ml 



Ml 



A, 



Ml 



ml Ami 



.1 TV 



Ml 



(9) 



(10) 



(11) 



The form factor v(k) and cr-mass M a — 1 GeV are fit to experimental phase shifts, as in Ref. [0. For free tt + -it~ 
scattering this force yields, when used in the T-matrix (see below), a scattering length which vanishes in the chiral 
limit, as it should. This feature induces off-shell repulsion below the 27r-threshold in spite of the fact that the positive 
Sq phase shifts indicate attraction. It is remarkable that the gap equation still shows a non-trivial solution, signaling 
pion pair formation, as we will show later. It is evident that bound pair formation, or even larger clusters of pions, 
can deeply influence the dynamics of the pion gas. 



In the sigma channel (£ = 0,1 = 0) we rewrite the gap equation ([[]) as 
A(fc) = -■ 



^3 (k,-k\ V i =q(E = 2(j) I k'-k') ^1 



(2tt) 3 N ' 1 ^ 1 ' ' 2E{k') 

With the form of our interaction solutions of this equation may be written as 



coth 



E(k') 
2T 



A(k) 



;(*) 



Si vi(k) + 5 2 v 2 (k) 



(12) 



(13) 



and Eq. ( |l2| ) reduces to two coupled non-linear equations for the "gap strengths" Si and 6 2 ■ For a non-trivial solution, 
one can show that <5i > —5 2 > 0. We also note that while A2 is always repulsive, \i(E) is attractive at E = 2/i only 
if I fi I > Ma-m^ /2y / 3Ml — 2ml ~ 40 MeV. This inequality is also the formal condition for a solution to the gap 
equation to exist at some temperature. Intuitively, we require at least some attraction because, as we shall see, a 
solution to the gap equation is connected to the existence of a pole in the T-matrix. We note that the repulsive part 
of the tt-tt interaction Eq. (^) helps to avoid collapse. This is different from our previous calculation jlij , which was 
performed with an entirely attractive interaction. The presence of this repulsion is a consequence of chiral symmetry 
and PCAC @- 

In the previous section, we introduced the critical temperatures T", at which the gap vanishes, and T^, where the gap 
has reached its maximum value and quasiparticle condensation occurs. Fig. 1 shows the numerical results for these 
temperatures in the fx-T plane. The T" (solid line) are obtained by linearizing Eq. (|j|): 
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Mk) = -\j^ (k,-k\V I=0 (E = 2^ | k',-k') coth^ , (14) 

while the Tjf (dashes) result when the gap strength increases to a point where E(k — 0) = 0, i.e., m, — // = <5i + 62- 
At high temperatures T > T" (region III), the system is in the normal state with no gap, while below the dashed line, 
T < (region I), there is macroscopic occupation of the k = mode. For Xjf < T < (region II), non-trivial gap 
solutions exist. Notice that for physically realistic solutions (T < 200 MeV, say) we have fi ~ m v , and w-m, <C m OT , 
and, in hindsight, are justified in neglecting relativistic corrections to the gap equation (see appendix). 

Fig. 2 shows the gap strengths Si (solid line) and —62 (dashes) versus temperature for a fixed chemical potential 
/i = 135 MeV. Again, we see that at high temperatures, in region III, only the trivial solution Si — 82 = exists. As 
the temperature drops to T" ^123 MeV, the order parameter A switches on, and we have a transition to a paired 
state in region II (see discussion below). Finally, at T = T~. ~ 77 MeV, the gap has reached its maximum value 
Si + S2 = m-K — /i ^ 3 MeV and quasiparticles condense in the lowest energy mode in region I. 

The change in the pion dispersion relation E(k) is investigated in Fig. 3 in the temperature range Tjf < T < T", 
for a fixed chemical potential of fi = 135 MeV. At T — T c u ~ 123 MeV (solid line), and above, we simply have 
the normal-state pion dispersion relation e(k) — u(k) — \i. With decreasing temperature the influence of the finite 
gap becomes visible at long wavelengths: The dot-dashed line shows E{k) for T = 115 MeV. A further drop in the 
temperature to T — Tjf ~ 77 MeV qualitatively changes the character of the pion dispersion relation to a linear, 
phonon-like dispersion at small fc. 



IV. THE IN-MEDIUM nn SCATTERING MATRIX 



We turn now to a discussion of the T-matrix Mi=q(E, K) for a pion pair with total momentum K =| K \ with respect 
to a thermal medium. Writing the on-shell T-matrix (c./. Eq. ( A.29| )) as 

2 

(k*\M I=0 (E,K)\k*) = J2 W*)vi(.k*)Ti(k*;s,K) , (15) 



where fc* 2 = s/4 — m 2 and s — E 2 — K 2 is the square of the total cm. energy, the Lippmann-Schwinger equation 



(A. 27) becomes a set of two linear equations for the functions 7^: 



2 

E 

3=1 



Sij - X j (s)g ij (s,K) 



i(k*;s,K) = Vi (k*) 



i= 1,2 



with 



9v(*>K) = \ j J^js < k ) TTTIa ^! + ta + , f 2 v 'J [l] 



(16) 



(17) 



ui(k) s — Alu 2 , + it] 

Here, (1 + /+ + /_) denotes an average over the angles of the cm. relative momentum of the pair. For thermal 



occupation numbers it is given by Eq. ( A.28|), and reduces to unity in free space and coth(o;(/c) — /i)/2T for vanishing 



total momentum K. We note that Eq. (|T^) does not incorporate the non-linear effect of the gap. 

The solid line in Fig. 4 shows | Aii = o | 2 for free space scattering. Compared to our previous calculation Q, the T- 
matrix is relatively flat above the resonance, this being due to the repulsion in the interaction at high energies. The 
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short dashes give | A4i = o | 2 in a thermal bath of T — 100 MeV and fi = 135 MeV, for K = 0. The medium strongly 
suppresses the cross section, an effect that also occurs in the (1=1,^=1) p-channel |@,||. At high cm. energies, the 
phase space occupation becomes negligible, and the cross section returns to its free space value. The three remaining 
curves show results in the same thermal bath, but for K — 200 MeV/c (long dashes), 1 GeV/c (dot-dashed) , and 
3 GeV/c (dotted). As K increases, the pair is boosted more and more out of the occupied phase space of the medium, 
and the cross section again returns to its free space value. We also see a threshold behavior in Fig. 4: as K becomes 
larger, a resonance peak emerges from below the threshold which continues to shift up in energy and strengthen until 
it coincides with the free scattering peak. We shall see below that this is the continuation of an upward shift of the 
Cooper pole in the T-matrix with decreasing phase space occupation |l4| ] . 

We consider now the existence of poles in the T-matrix, first for the special case of zero total momentum, K = [[lJJ , 
and define the determinant function 



F„, T (E) = -det 



5 %3 - X J (E)g l] (E) 



(18) 



This function is shown in Fig. 5 for five different temperatures (solid lines) at a fixed pion chemical potential of 
135 MeV. The intersection of these curves with zero (horizontal dashes) below 2m w (the bound state domain) 
gives the pole position. We see that a pole always occurs provided the temperature lies above some critical value 
Tq w 47 MeV, for which the pole is at threshold and ceases to exist. This Tq is close to the lower critical temperature 
for the gap, T^, where the excitation spectrum vanishes at k = and quasiparticles begin to condense as singles. Thus, 
the bound state and gap solution disappear at a similar critical temperature; differences are ascribable to the fact 
that we use free quasiparticle energies in the T-matrix. 

There is a second special temperature Tq, for which a pole exists at E = 2fi (see Fig. 5). It is identical to the upper 
critical temperature T" at which the gap vanishes, as may easily be seen by rewriting the T-matrix for E near 2\x, 

(k,-k\T I=0 (E)\k',-k') = Z{k)—^—Z{k') . (19) 

hi — ZfJ, 

In the non-relativistic limit, Z(k) follows as (see appendix) 

Z( k ) = -- j fX^{k,-k\V I - Q (E = 2 l i)\k l ,-k l ) , y N coth Cj(fc/) ~ M , (20) 
y ' 2 J (2tt) 3 V ' 1 >y 1 2(uj{k r ) - fi) 2T " ' V ' 



which is precisely the same condition as for T", Eq. (|14[). The gap equation (|12|) thus reduces to the T-matrix pole 
condition at the particular energy E — 2/i. In fermion systems, this is the well-known Thouless criterion |l6| for the 
onset of a phase transition to a pair condensate. We note that the Thouless criterion is only approximately valid if 
relativistic corrections are included. 

Several observations can be made. Firstly, one always obtains a pole in the T-matrix if the temperature lies above 
Tq(h). Thus, at fixed [i, no matter how weak the interaction strength is (provided it is attractive in the neighborhood 
of the 2tt threshold), one always obtains a pole for sufficiently high temperatures (for fermions at a sufficiently low 
temperature). In practice, Tq (and Tq) will exceed sensible values for pions as soon as fj, drops below ~ 130 MeV, 
since they are increasing functions of fi. Secondly, for a fixed interaction strength, the pole position shifts downward 
with increasing temperature (for fermions: pole position moves up with increasing temperature). As a function of 
temperature, we therefore see a behavior for bosons opposite to that for fermions. 
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The fact that increasing temperature reinforces the binding is somewhat counterintuitive, but is an immediate con- 



sequence of the coth-factor associated with bose statistics in Eq. ( A. 27 ). Indeed, one realizes that the coth-factor 
increases with increasing temperature and thus effectively enhances the two-body interaction. We can therefore always 
find a bound state for arbitrarily small attraction: it suffices to increase the temperature or, equivalently, the density 
accordingly. This is opposite to the fermion case where the corresponding tanh-factor suppresses the interaction with 
increasing temperature. Therefore, in the fermion case, even at the T-matrix level there exists a critical temperature 
where the Cooper pole ceases to exist. For bosons, on the other hand, once one has reached T = Tq, a bound state 
(here E^ < 2m n ) exists and the bound state energy simply continues to decrease as the temperature increases. Of 
course, this becomes unphysical as soon as the density of pairs becomes so strong that the bound states start to 
obstruct each other, and finally dissolve at an upper critical temperature (Mott effect). Precisely this non-linear 
effect is very efficiently taken care of in the gap equation. In spite of the fact that we still have a coth-factor in the 
gap equation, there is now a crucial difference: the argument of the coth-factor is the quasiparticle energy, Eq. (Q), 
(over T) and thus, due to the presence of — A 2 (k) in E(k), the origin of the coth is shifted to the right with respect to 
the T-matrix case. Now, as T increases, the only way to keep the equality of the gap equation is for A(fc) to decrease 
- this pushes the origin of the coth back to the left, counterbalancing its increase due to the increasing temperature. 
Of course, this only works until A = 0, i.e., until the temperature has reached T c u . This is precisely the temperature 
Tq for which the bound state in the T-matrix reaches an energy E 2 ^ = 2(m w — We therefore see that in spite 
of the fact that the bosons prefer high phase space density, the formation of bound states ceases to exist beyond a 
critical temperature - just as for fermions. 

Lastly, we return to the behavior of the pole for varying total momentum K , and the threshold effect seen in Fig. 4. 
Since F^^t(s,K) becomes complex above threshold, we show in Fig. 6 its magnitude for fixed T = 100 MeV and 
/i = 135 MeV, and various values of K. As expected, for increasing K (i.e., decreasing phase space occupation felt 
by the two pions in question) the pole (zero of F) moves up in energy until it dissapears at some critical momentum 
100 McV/c < K c < 250 MeV/c. For K > K c , the now non-zero minimum of the determinant function continues to 
shift to higher energies, corresponding roughly to the similar shift in the resonance peak in Fig. 4. 



V. DISCUSSION AND CONCLUSIONS 



In the previous section, we investigated the effect of a thermal medium on the pion dispersion relation at low momenta 
k. In particular, one finds a critical temperature Tjf at which the pion dispersion relation is linear (phonon-like) in 
k for small k. This result is independent of the details of the interaction and characteristic of any bose system (see 
Ref. [0 for the case of 4 He). 

Such a change in the pion dispersion relation at low temperatures would influence the pion spectra at low momentum. 
For this to occur, rather large medium phase space occupation numbers are required. In particular, for a physically 
reasonable system with, say, T < 200 MeV, this means that we require large chemical potentials. In fact, dynamical 
calculations show that a buildup of /i can indeed occur, provided that the scattering rate is sufficiently large 
compared to expansion rate and the inelastic collisions have ceased to be a factor. 

To demonstrate the possible effect in a qualitative way, consider the pion transverse momentum spectrum for longi- 
tudinally boost invariant expansion 
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ssss " MV) <#^ «*?>*.<tt> • (21) 

where _ftTi is a Mc-Donald function, mt is the transverse mass, and the normalization volume itR 2 t is of the order 



200 to 300 fm 3 |18|. At mid-rapidity, the transverse mass coincides with the full energy of the pion, and we follow 
Rcf. 1 19 1 in replacing mt by the in-medium pion dispersion relation E(k) + \i derived in the previous section. Of 
course, as remarked in Ref. p9| , the use of this procedure is rather tenuous since the system is by definition still 
far from freeze-out. In a dynamical calculation, hard collisions would re-thcrmalize the system at ever decreasing 
temperatures. 

In Fig. 7, the thermal transverse momentum spectrum for pions with /j = 135 MeV and T — 100 MeV > T~([i) 
is shown with (solid line) and without (dashes) the effect of the gap energy. Essentially, the presence of a large 
chemical potential gives the spectrum the appearance of one for small-mass particles, and the gap energy, which 
causes E(k) ~ k for long wavelengths, strengthens this effect. We would like to mention here again that the use of 
our force, Eq. (||), which respects chiral symmetry constraints [jl7|, considerably reduces the effect of binding with 
respect to a purely phenomenological interaction, fitted to the phase shifts (see, for example, Ref. 0]). This stems 
from the fact that the expression (||) becomes repulsive sufficiently below the 2m w threshold. This is not the case for 
commonly employed phenomenological forces fl7j| . As a consequence, the effect we see in Fig. 7 is relatively weak, but 
one should remember that the force (j^) is by no means a definitive expression. It is well known that in a many body 
system screening effects should be taken into account. Whereas in a fcrmi system this tends to weaken the force, it is 
likely that screening strengthens it in bose systems. In this sense our investigation can only be considered schematic. 
A quantitative answer to the question of bound state formation in a hot pion gas is certainly very difficult to give. 
Qualitatively, the curves in Fig. 7 agree with the trend in the pion data at SPS to be "concave-up," but this is 
mainly an effect from the finite value of the chemical potential ||,[l8| . While the gap changes the spectrum by a factor 
of ~ 3 at m t — run ~ 0, this region is not part of detector acceptances. 

In summary, we have shown that finite temperature induces real poles in the 2ir propagator below the 2m w threshold, 
even for situations where there is no 27r bound state in free space Q . The situation is analogous to the Cooper pole 
of fcrmion systems, and we therefore studied the corresponding bosonic "gap" equation. This equation has non-trivial 
solutions in a certain domain of the \i-T plane. Such a region always exists, even in the limit of infinitesimally weak 
attraction. This is different from the T — case discussed by Nozieres and Saint James jl2j , where a nontrivial solution 
to the gap equation only exists when there is a two boson bound state in free space. Our study has to be considered 
preliminary The final aim will be to obtain an equation of state for a hot pion gas within a Bruckner— Hartree- 
Fock-Bogoliubov approach. Also, the subtle question of single boson versus pair condensation must be addressed (see 
Ref. Jl^ ] and references therein). Furthermore, the fact that we obtain two pion bound states in a pionic environment 
leads to the speculation that higher clusters, such as four-pion bound states, etc., may also occur, and perhaps even 
be favored over pair states. Such considerations, though interesting, are very difficult to treat on a quantitative basis. 
However, substantial progress towards the solution of four body equations has recently been made and one may 
hope that investigations for this case will be possible in the near future. 

We are grateful to P. Daniclcwicz for useful discussions. This work was supported in part by the U.S. Department of 
Energy under Grant No. DE-FG02-93ER40713. 
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VI. APPENDIX: DERIVATION OF T— MATRIX AND GAP EQUATION 



This appendix is devoted to a derivation of the gap equation for a bosonic system governed by a field-theoretic 
Hamiltonian. The basic problem one has to deal with is the formal introduction of a chemical potential for bosons, 
since the total bosons number operator (i.e., n + + 7r~ + 7r°) does not commute with the Hamiltonian. Hence, if we 
consider a pion gas at a typical temperature of 200 MeV, it will correspond to zero chemical potential. However, for 
a system lifetime on the order of tens of fermi, the inelastic collision rate is negligible. Therefore, provided the elastic 
collision rate is sufficiently large, a thermal equilibrium with a finite chemical potential may well be reached. 

Let us consider a pion system at temperature T. Inspired by the linear cr-model, with form factors fitted to the tt-tt 
phase shifts, we take the Hamiltonian 



H = Hn 



H 



int. 



(A.l) 



where Hq is the kinetic Hamiltonian for the ir and "cr" mesons 



(A.2) 



The index "1" refers to the momentum and isospin of the pion, and "a" to the momentum and identity of the heavy 
meson carrying the interaction. The interaction Hamiltonian has the form 



(a a +al a ) (b\b\ + b^ 2 ) 



(12 I W I a) 



b\blb 3 bi + -{b- 1 b^ 2 b 3 b 4 + blblbl 3 bt i ) 



(12 | V | 34) 



(A.3) 



In the linear cr-model one has (L 3 is a normalization volume) 



(12 | W | a) 
(12 | V | 34) 



1 



1/2 



2wiL 3 2uj 2 L 3 2fl a L 3 
1 

2wiL 3 2uo 2 L 3 2uj 3 L 3 2ui A L 3 
Ml - m 2 r 



v(k* 12 ) (2TT) 3 8(k l + k 2 
1 1/2 



K n ) : di2 



v(k\ 2 ) v(k* 34 ) (27r) d <J(fei + k 2 



ki) 



■I TT 



2ujiuj 3 ~ 2ki ■ k 3 ~ ml 2wi^ 4 - 2fci ■ fc 4 - ml 1 

012 034 + 013 024 77^ h 0l4 023 



Ml 



Ml 



(A.4) 



(A.5) 



The form factor taken at cm. momentum k* of the pion pair is fitted to the experimental phase shifts, and M a is 
the cr-mass. The static quartic interaction contains the tt 2 tt 2 interaction of the cr-model, and the t and u channel 
tr-exchange terms. We neglect the t and u dependence of the denominator (see Ref. fT7| ), since their effect is extremely 
small. Further, terms like ab'b and b^bbb have been dropped, since they are not essential for our purpose. 

The Dyson equation for the pion propagator 

We now derive the equation of motion for the pion propagator 

Gil(M') = ( -iT(h(t)b\(t') 



(A.6) 



where the b\ are normal Heisenberg operators bi(t) = exp(iHt)bi(0) exp(— iHt). In principle, the extension to finite 
temperature requires a matrix formulation (real time formulation) or Matsubara Green's function. However, for 
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simplicity we consider here the normal zero temperature G, and replace it by a thermal propagator at the end. We 
have checked that the final result is not modified. 



After standard manipulation, using the Hamiltonian ( A.l ), we obtain 

(i^-c^Gnftf) = 5{t-t') + J dt"Xi(t,t")Gn(t",t') 

with £i(f, t') — Sf (t, t') + t'). The static part of the mass operator is 

E?(t,f) = Yl (4&2> (12 I V | 12) - , 



(A.7) 



(A. 



T([ffi nt) 6i](t) [ff ta t,6 1 ] t (< / ; 



(A.9) 



while the dynamical part is given by 

Ef(t,0 = 

Making a standard factorization approximation, we obtain 
£?(M') = *E G 2 2 (M') 

2 

x ^ - iT^ (12 | W | a) (a a + al a )(t) + ^(12 | V | 34) (b 3 b 4 + bl 3 b f _ 4 )(t) 
^l+<rl a ,)(t') (a'\W\ 12) + \{b\,b\, + &_ 3 ,&_ 4 ,)(f) (3'4' | V | 12)])^ 

with G2 2 (M') = (-iT(bl(t),b 2 (t'))), and there is an implicit summation over repeated indices. 
Extraction of the condensates 

In the above expression the operator 6364 connects states with N particles to states with N + 2 particles. Among 
these states those with excitation energy 2/j, play a prominent role (Cooper poles). To separate the influence of these 
states, we split the fluctuating part of the operator from the condensate 



(A.10) 



b\b\{t) = (&T&T(t)) + : b\b\{t) 



(All) 



The time evolution is 



{b\b\(t)) = (fttftL,) e ^ <5 3 ,_ 4 , 



where (&3&I3) is the usual time independent pion density. Similarly, we obtain 

b 3 bi(t) = (636-3) e"^ 4 <5 3 ,_ 4 + : b 3 b 4 (t) : . 

We now extract the condensate part of the cr-field operator from the fluctuating part: 

a a (t) = (<r a (t)) + s a (t) . 

The equation of motion gives 



(A.12) 



(A.13) 



(A. 14) 
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i^{cr a {i)) = tl a (<r a (t)) + \{{bib 2 + &L1&L3) (a I W | 12) 



We look for a solution of the form 



(a a (t)) = (ie- 2 "' + Be i2 "V«o 



(A.15) 



(A.16) 



A and £? are straightforwardly obtained from the equation of motion: 



M«)) 



1 (hb-jjoiwii-i) e _ 42Mt _ 1 (^(Ol^ll-l) ^ t 

2 M CT - 2(i 1 2 M ff + 2/i 



(A.17) 



In the expression of the dynamical mass operator one can extract a Cooper pole part, where only the condensates 
occur. The remaining part involves only the fluctuating pieces. Grouping the latter with the static mass operator, 
we can write 



Ei(M') = Sic(M') + Sm(M') 



(A.18) 



where £i#(i,i') is the normal Hartree mass operator which depends on the full in-medium T-matrix: 

Eih(M') = G 22 (t,t') (12 | T(M') I 12) ■ 



(A.19) 



with 



(12 | T(t,t') | 34) = (12 | V | 34) 5(t - t') 
■iT 



(12 | W | a) (s a + sl a )(t) + | (12 | V | 56) (: 6 5 6 6 + &1 B &! 8 :)(i) 
(4+«-a')(*0 («' I W | 34) + i(: bl4, +&_ 5 '&-6' 0(f) (5'6' | V | 34)])^ 



(A.20) 



Using the Dyson equation for the b and s operators, it is a purely technical matter to show that this scattering 
amplitude satisfies a Lippmann-Schwinger equation. In energy space, and in the 1 = channel, it reads: 



(12 | T J=0 (E) | 34) = (12 | V I=0 {E) | 34) + 1 (12 | V I=0 (E) | 56) G 5 2 l(E) (56 | T I=0 {E) | 34) 



where G 2 - K {E) is the in-medium 27r propagator 
Gill. El ' 



1 



E — (u 5 + cj 6 ) + ??7 E + (cJ5 + oj 6 ) + Z77 



(1 + h + h) , 



(A.21) 



(A.22) 



with thermal occupation numbers f(k) = [exp(ui(k) — (i)/T — 1] . As mentioned above, we have checked that using 
the correct matrix form of the two pion propagators instead of (A.22) yields the same final result. In Eq. ( A.2l| ), 
Vi = o(E) is the effective tt-tt potential in the 1 = channel which incorporates all the tree level diagrams. For total 
incoming momentum K = ki + k 2 = k% + k^, it reads 

. . / 1 X 1 / 2 - . 

(fci,fc 2 | V I=0 {E) | fc 3 ,fc 4 ) = - - - - (fci,fc 2 |M B (£)|fc 3 ,fc4> , (A.23) 



2wi2tj22cL;32ci;4 



where the bare invariant interaction Mb is 



12 



(ki,k 2 | M B {E) I k 3 ,k 4 ) = (k* 12 | M B (s) I fc 3 * 4 ) = ^(s)vi(k* 12 )vi(k 



34) j 



(A.24) 



with 



Vi(k) = v(k) = [1 + {k/8m„f 



V2{k) = ^l v{ k) 



X 1 (s) 



n 

■I 7T 



s - mi 2mi 



Ml Ml 



_ Ml- ml Ami 
2 " ft W 



(A.25) 
(A.26) 



In these equations s = E 2 — K 2 is the square of the total cm. energy, and the k*j are the magnitudes of the relative 
3-momenta in the cm. frame 



* 2 2 , j * 2 ± 

= m^ + k i3 = - 



> i,i = 1,2 or 3,4. 



The form factor v(k), Eq. (A.25), and c-mass M a = 1 GeV have been fitted to the experimental phase shifts. 



The Lippmann-Schwingcr equation for the invariant T-matrix, Aii=o, may finally be rewritten in a form suitable for 
practical purposes: 



(k* u | M I=0 (E,K) | fc* 4 ) = (k* 12 | M B (s) I k* 34 ) 



+ 



1 fd 3 k 



(k* 12 | M B (s) I k* 6 ) 



1 <! + /+ + /- 



^{^I^^)!^) • (A.27) 



2 J (2tt)3 — 1 ' — ' 1 u 5 * 6 fl -4a& a + w? 
In the special case of a single fireball of temperature T and chemical potential /i, the angle average factor is given by 
(1 + U + /-> = in ^{1° + - , (A.28) 



7/3* s - 8 sinhlhK, -/)!%*,) -/i}/2T] ' 

where /3 and 7 are the velocity and gamma-factor of the pair with respect to the bath. This factor reduces to 
coth(wg 6 — fJ.)/2T for vanishing incoming total momentum K. 



Eq. ( |A.27| ) is solved by a separable ansatz 

2 

(k* 12 I M X=0 (E,K) I fc 3 * 4 ) = J2 X *( s ) v l (k* 12 )r l (k* 34 ; S ,K) , 

i=l 

where the functions r, obey the coupled set of equations 

2 

Tj(k;s,K) = Vi(k) , i=l,2 



(A.29) 



E 



(A.30) 



with 



'—^3 "tix ~ + { + .2 + , 4r ^( fc ) 



(27T) 



(A.31) 



The gap equation 

To obtain the Cooper piece of the mass operator we must simply replace a and bb by (a) and (bb) . According to the 
previous result, we find after some straightforward algebra, and noting that the index 2 is necessarily —1, 
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The important point is that F± involves the I = 0, £ = energy dependent 7r-7r potential at E = 2/i: 

1 /" d 3 fc 2 



| V} =0 (F = 2/i) | fc 2 ,-fc 2 ) (626-2) 



2 J (2tt) 3 

In energy space, Eic(w) is 

S lc (o;) = F 1 2 Jdre™{9{T) (bt(r),h (0)) + 0(-r) <M0), b|(r))} 
x ( e -^ r + c l2 ^ T ) . 

Taking for 6{(r) the bare time evolution 6{(r) = 61 e 1 " 1 ", and keeping only the real part, we finally find 

1 1 



(A.32) 



(A.33) 



\ui + u>i — 2/i £J + wi + 2// 
The first term is the usual non-relativistic result, and the second one corresponds to a relativistic correction. 



(A.34) 



(A.35) 



Reinserting the result (A.35) into the Dyson equation for the pion propagator, and ignoring the Hartree correction, 
we find that the pole of the pion propagator is the solution of 

(oj - fl) + [uj\ - fl) 



(u - fif = (wi - m) 2 - Ft 



1 



(uj — fl) + (uii - fl) + ifl 



(A.36) 



The second term in the square brackets represents a relativistic correction to the standard dispersion relation, since 
for typical non-relativistic situation one has 



[i ~ m n , u> — TOtt <C m„, and u>\ — <C m v 



(A.37) 



Calling 



A 2 = Fl 



1 



Ei + (wi - fl) 



Ei + (lui - fi) + ifi 

the quartic equation can be approximated by a quadratic one in terms of ui — fi: 



E 



with a gap equation following from Eq. ( A.33| ) and ( A.35Q 

Ai = -- [ (kx,-ki I Vi- (E = 2u) I k 2 ,-k 2 ) — coth — 



Fi + (loi - fl) 



Ei + (uji - fi)+Afi 



1/2 



(A.38) 



(A.39) 



(A.40) 



which is the standard gap equation with a relativistic correction. The presence of the factor 1/2 is somewhat uncon- 
ventional, but is simply related to the fact that the matrix element of the interaction incorporates the exchange term. 
Note that the factor 1/4 in front of the quartic term of the interaction Hamiltonian Eq. ( |A.3| ) has the same origin. 

To calculate the occupation number, we note, using the explicit form of that the bare vacuum is the vacuum of 
quasi-particle operators Bi , such that 



14 



2Ei 



1/2 



2Ei 



1/2 



(A.41) 



Using (B^B) = [exp(Ei/T) - l]" 1 , it follows that 

(b\bi) = 



2E 1 



^coth^ 
2T 



(6i6_i) = — - coth^- 
x ' 2E X 2T 



(A.42) 



The Thouless Criterion 

We may also obtain the condition for having a pole in the T-matrix at E — 2/i. For E near 2/i, we write 

1 



(ki,-k! | T T=0 (E) | k 2 ,-k 2 ) = Z 1 



E-2n 



Multiplying (A. 27) by E — 2/i and taking the limit £7 — > 2/i, one obtains an equation for Zi 

w 2 - /i 



Zi 



— ^ (fei,-fc! | V /= o(£ - 2/i) | fc 2 , -fe,) 2 



coth ■ 



2T 



1 



£ 2 + (u>2 - /x) 
£" 2 + (w 2 - /i) + 4/i 



(A.43) 



(A.44) 



In the non-relativistic limit, i.e., neglecting the last term in the square brackets, this equation coincides with the 
linearized form of the (non-relativistic) gap equation. This is just the Thouless criterion: the gap equation begins to 
exhibit non-trivial solutions at the point where the T-matrix has a pole at zero energy, (H — fiN) = E — 2/i = 0. We 
see that the Thouless criterion is only approximately valid if relativistic corrections are included. 
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FIGURE CAPTIONS 

Fig. 1 The critical temperatures (solid line) and (dashes). 

Fig. 2 The gap strengths 8\ (solid) and — &2 (dashes) vs temperature T , for fi = 135 MeV. See text. 

Fig. 3 The quasiparticle dispersion relation at fixed fj, = 135 MeV, for three temperatures: T~ (dashes), T = 115 MeV 
(dot-dash), and T c u (solid). 

Fig. 4 The square of the on-shell invariant T-matrix for / = 0,£ = in free space (solid line). Also shown are the 
results in a thermal bath with T = 100 MeV, // = 135 MeV for total momentum K = (short dashes), 200 MeV/c 
(long dashes), 1 GeV/c (dot-dashes), and 3 GeV/c (dotted). 

Fig. 5 The T-matrix pole function T Mj t(£0 at various temperatures for /i = 135 MeV and K = 0. Here, Tq = T" « 
123 MeV, T l c « 77 MeV, and T £ w 47 MeV. 

Fig. 6 The magnitude of F^ T (s,K), Eq. |[|), for a thermal medium with T = 100 MeV and fx = 135 MeV. The 
solid, short dashed, long dashed, and dot-dashed lines correspond to K = 0, 100 MeV/c, 250 MeV/c, and 500 MeV/c, 
respectively. 

Fig. 7 The thermal pion transverse mass spectrum at midrapidity, for T = 100 MeV and fi = 135 MeV, without the 
effect of a gap (dashes) and including the effect of a gap (solid line). 
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